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^ . 1. Introduction 



Following up the investigation of holonomy in the Schwarzschild spacetime by 
Rothman et al |jl|, Maartens et al |^ have extended the work in ||l| on closed 
circular trajectories to the Kerr spacetime where the spacetime rotation further 
complicates matters by introducing new features. Maartens et al ||^ study the parallel 
transport of a generic vector X {u in their notation) in the equatorial plane of the 
Kerr spacetime along the orbit of the spacelike Killing vector associated with the 
axial symmetry. Assuming the standard Boyer-Lindquist form of the Kerr metric 
with (x'^ , , x'^ , x^) ~ {t,r,6,(j)), the time coordinate t parametrizes each orbit of 
the stationary symmetry Killing vector field dt and the azimuthal coordinate </> 
parametrizes each orbit of the axial Killing vector d^. The transport equations reduce 
to a single fundamental second-order ordinary differential equation for the 0-derivative 
of the ^-component of X 

at/) aq) r r 
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Recalling that r is a constant along these orbits, this is trivially solved for Z, resulting 
in solutions characterized by a real or purely imaginary dimensionless frequency /. 
The other components , X'^ and X^ , and X* itself are then easily calculated. 

Maartens et al |^ discuss the holonomy question by studying the behavior of 
these solutions, which clearly depends on the sign of the frequency-squared function 
= F/r^ , with the roots of the quartic equation F = in r distinguishing oscillatory 
and exponential solutions of the system. They show that F has exactly one positive 
root r* (called rcrit in the Schwarzschild discussion of Rothman et al but find 
no obvious significance for this value. When F > 0, which occurs for r > r, , 
rational values of the real frequency / lead to the transported vector X returning 
to its original value after a certain number of loops. Here we show that and its 
sign have a geometrical significance directly associated with clock effects |^ and then 
extend the study to the closed geodesic circuits associated with these clock effects and 
also studied by Rothman et al for the Schwarzschild spacetime. The results hold 
for any stationary axisymmetric spacetime for which clock effects can be studied. This 
requires the existence of pairs of oppositely-rotating timelike circular geodesies which 
intersect periodically, forming a double helix in spacetime contained within a single 
t-(f) symmetry group orbit. 

For a general family of stationary axisymmetric circularly-rotating observers, with 
an intermediate angular velocity between the two distinct angular velocities of such 
a pair of timelike circular geodesies within a given circular orbit cylinder, if the two 
oppositely-rotating geodesies (as seen by these observers) depart from a given observer 
world line simultaneously, then they will in general not return to that world line 
simultaneously. Instead there will be a proper time difference in their arrival times 
as measured by the observer clock. This is the observer-dependent single-clock clock 
effect, where the observer both defines the complete closed circuit and measures the 
difference in arrival times for it, interpreted as the difference in the observer- measured 
orbital periods. The observer-dependent two-clock clock effect is instead the difference 
in proper times elapsed on the clocks carried by the geodesies themselves for the same 
observer-defined circuit. Finally the observer- independent two-clock clock effect is 
the difference in proper times elapsed on the clocks carried by these geodesies between 
three successive meeting points; this effect is equivalent to the observer-dependent two- 
clock clock effect for the 'geodesic meeting point observers' (GMPOs) when the latter 
are timelike. The trajectories with constant angular velocity which pass through every 
other meeting point of the pair of geodesies are the geodesic meeting point (GMP) 
trajectories. When these trajectories are timelike, they define the world lines of the 
GMPOs |,|. 

It turns out that the zero of the frequency-squared function f^ corresponds to 
the radius ^(gmp) of the observer horizon for these GMPOs. Each such observer sees 
the pair of corotating and counter-rotating geodesies which leave it simultaneously 
also return simultaneously, and thus sees vanishing single-clock clock effect. This 
observer family has a number of special geometrical and kinematical properties and 
is also useful (by definition) in the discussion of the observer-independent two-clock 
clock effect. Outside the radius r*, these observers are able to 'resist' the clock effect 
by counter-rotating sufficiently fast to synchronize the geodesic arrival times, but 
within this radius the rotation of the spacetime itself is too great to overcome. The 
GMP trajectories are well-defined whether the circular geodesies are timelike, null or 
spacelike; it is convenient to define their direction using their unit tangent when the 
tangent is non-null. Only when timelike can these GMP trajectories be interpreted as 
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the congruence of world lines of a test observer family. 

The clock effect compares the spacetime arclength for different pairs of timelike 
paths between two spacetime points, each pair together forming a closed loop in 
spacetime. One can also compare parallel transport around such closed loops, 
investigating special elements of the holonomy group of the spacetime. This requires 
first solving the equations of parallel transport for a vector along a general circular 
orbit with constant angular velocity, which is done here. For comparison, the simpler 
case of closed (/)-coordinate curves is studied, followed by the loops formed by the two 
oppositely-rotating timelike geodesies. In the first case the geometry depends crucially 
on the GMPOs, while in the latter case both the GMPOs and the Carter observers 
come into play. The geometry of circular orbits and of the circular geodesies which 
give rise to both of these observer families requires preliminary study. The extremely 
accelerated observers [|[ ^ also appear in this discussion as the only observers 
which see the parallel transport geometry symmetrically for pairs of orbits with 
equal magnitude but oppositely-signed relative velocities like the circular geodesies 
themselves. 

The holonomy invariance found at a dense subset of special radii and particular 
loop numbers for (/)-coordinate loops and geodesic loops in the static case [0 and for 
0-coordinate loops in the stationary case and referred to as a 'band' structure 
is reduced to a set of radii at most of measure zero for closed geodesic loops in 
the nonstatic case. In the Kerr spacetime the dimensionless frequencies associated 
with parallel transport rotations can be expressed as ratios of the proper and average 
coordinate periods of the circular geodesies. Establishing these results is an elegant 
exercise in geometry. The process of discovering them, however, initially required 
a computer algebra system to jumpstart the investigation and to derive certain 
consequences which could then be confirmed by hand. The entire discussion is a further 
example of how various geometrically defined circularly orbiting observer families 
continue to play a key role in understanding the geometry of stationary axisymmetric 
spacetimes. 

2. Metric and circular orbits 

Consider the spacetime metric in symmetry-adapted coordinates 

ds^ = ds(t_^) -I- grrdr^ + geedO^ , (2.1) 

with the metric of the circular orbit cylinder written respectively in coordinate and 
lapse-shift notation 

ds(t,0) = gttdt'^ + 2<?t</.dM0 -t- 500d(/)2 = -NMt^ + 5^^(d(/) + N'^'dtf , (2.2) 

where in the case of Kerr spacetime, {t, r, 9, (j)} are the usual Boyer-Lindquist 
coordinates adapted to the family of static observers, whose world lines are along the 
time coordinate lines, and to the zero angular momentum observers (ZAMOs), with 
4- velocity n = N~^{dt — N'^'d^) and world lines orthogonal to the time coordinate 

— 1/2 

hypersurfaces. Let = g^^ di be the orthonormal frame vectors along the time 
slices (to be referred to as the 'spherical frame'), where i = r,9,(j). dt and are 
Killing vector fields associated respectively with stationarity and axisymmetry; 
has closed orbits cj) : 2tt. All considerations below will be limited to the region 
of spacetime where the 1-form dt remains timelikc and this lapse-shift decomposition 
remains valid, namely outside the (outer) horizon r = ?'-|_ in the Kerr spacetime, the 
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(outer) surface where N — (— 0. Static spacetimes occur as the special case 
gt4> = = N'^. Exphcit values of the various metric quantities for the Kerr spacetime 
are given in Table 1 of 

1/2 

The quantity R — g^^ represents the circumferential radius of the (/i-coordinate 
circles, converting to arclength so that the total circumference of such a circle is 
2-kR. On the other hand 

n{4>^ nf = -9~'/'Mg'/;)ir = -(Ini?),, = -1/p (2.3) 

is the (signed) Lie relative curvature of the 0-coordinate lines with respect to the 
ZAMOs [y and p is the corresponding radius of curvature, assuming that the radial 
curvature component is negative, which is true outside the horizon in Kerr. The ratio 
R/p = R^f is 1 in the flat spacetime case and the r-(f) surface is just a flat plane in 
Euclidean TZ^. When this ratio satisfies < R/ p < 1, the r-(/) surface can still be 
embedded in Euclidean TZ^ ||, |^. The tangent cone to a given (/)-coordinate circle 
then has a positive deficit angle 

A(def) = 2tt[1 - R/p] ^R/p=l- A(dcf)/(27r) (2.4) 

which helps gives a physical interpretation to the rotation induced by parallel 
transport after one circuit around the ^-coordinate circle as described below. For 
the Kerr spacetime, one has < < 1 outside the horizon, with R/p increasing 
monotonically from its limiting value at the horizon where g^^ ^ to its asymptotic 
value 1 as r oo ||]. 

The trajectory of a Killing vector field r] — dt + (d^ is a circular orbit whose 
tangent vector, when non-null, can be normalized 

[/(C) = r(C) [dt + Cd^] = 7(C) [n + lyiOe^] , (2.5) 

where i^(C) is the ZAMO relative velocity in the direction and the ZAMO gamma 
factor is 

liC) = \1 - HOT'^' = NT{C) , (2.6) 
while the coordinate gamma factor r(C) is defined by 

- (dt + Cd^) ■ {dt + Cd^,) = - - 2g^tC - 9tt = -g^^C - C-)(C - C+) 

^ e(C)'r(C)-2 , (2.7) 

and e(C) = [—U{C)-U{Q]^^'^ is respectively 1, 0, i for timelike, null and spacelike orbits. 
The zeros C± of this function are the angular velocities of the null circular orbits; when 
C = C±, then e(C) = and the rescaling factor T{(±) is arbitrary. When C- < C < C+ 
so that e(C) — 1, the quantity r(C)~"'^ converts intervals of coordinate time into proper 
time elapsed along these timelike trajectories 

At ^TiCy^At^jiCy^NAt . (2.8) 

In the non-null case one can also introduce the orthogonal unit vector with angular 
velocity C, defining the bar map 

UiO = UiC) = 7(C) [KC)« + e^] , U{0 ■ UiO = eiCr , 

DiO = KC) = 1/KC) , (2.9) 

which is the reflection across the null direction with the same sense of rotation as U{C) 
with respect to the ZAMOs. The corresponding one-form is 

UiCf = 7iOY{Cf , YiC? = - Cdi) . (2.10) 
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where b indicates index lowering. When U{() is timehke, U{C) is the unit vector in 

the (/(-direction in its local rest space and {e^, e^, U{C)} is the boosted spherical frame. 

The relationship between the angular velocity ( characterizing a circular orbit 
and the ZAMO relative velocity is 

Z/(C)=iV-^i?(C-C(„mp)) ^ /^(nmp) =KC(nmp)) =0 , (2.11) 

where 

C(nmp) = i(C+ + C-) = -g<^t/g^4> = -N'^ (2-12) 

is the angular velocity of the ZAMOs ("null meeting point observers"), whose world 
lines contain the meeting points of the oppositely rotating null orbits. 



3. Circular geodesies 

The circular geodesic 4- velocities 

f^± = r±(9t + C(gco)±50) = 7±(" + i^±e^) (3.1) 

are parametrized by the angular velocities C(geo)± = U^/Uj., which are determined by 
the radial component of the geodesic equation 

-^+T''^pU^Ui=0 . (3.2) 
Since U±. = dx" /dr, the radial component is 

^+T-c.^U^Ui=0 , (3.3) 
and setting r equal to a constant gives a quadratic condition for the angular velocity 

I"'00Cfgeo)± + 2r'',^tC(geo)± + r'^tt = (3.4) 

equivalent to 

500,rC(geo)± + 2.90(,rC(gco)± + gtt,r = . (3.5) 

It is assumed that these two roots are distinct and satisfy C(geo)- < < C(geo)+, which 
is the case in the Godel spacetime in its usual stationary axisymmetric form but not 
in Minkowski spacetime expressed in rotating coordinates, for example. 

In order to have circular geodesies at all, it must be true that the second 
nonignorable coordinate can also be constant along them, leading to a condition on 6 
once it is set equal to a constant ^(o) in the corresponding component of the geodesic 
equation 

0+r%C/2C/|=O. (3.6) 

This condition 

94>4>,oC{geo)± + '^94>t,eC{seo)± + 9tt,e = (3.7) 
requires that these 6 derivatives of these three metric components vanish at some 
value of 0. This is true for the equatorial plane = 7r/2 in the Kerr spacetime 
because of its reflection symmetry — > 7r/2 — ^, or for any constant value of ^ in a 
cylindrically symmetric spacetime if do is chosen to be the Killing vector associated 
with the additional translational symmetry, like the Godel spacetime, for example. 
It is also true for the equatorial plane of any stationary axisymmetric spacetime 
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C(ginp) - (C(goo)+ +C(gco)-)/2 


27r/T'(o) - (C(goo)+ - C(gco)-)/2 
= '^(0)7(1 - a^^^o)) 


^(car) = (^(gco)+ +^(goo)-)/2 

-a - (T+-T_)/2/(27r) 


T^(0)/(2^) = (C(-L)+-CteL)J/2 



Table 1. Table of angular velocity and period parameters associated with the 
circular geodesies. The Carter observers are associated with the static observer 
clock effect. 



that satisfies the condition of reflection symmetry about the equatorial plane, i.e. 

^5 = 9ti}^{t^ r, TT — 0, 0); in fact, the partial dervitative of this condition with 
respect to 6 implies that g^j^^e — at 6 ^ tt/2. It is not true in the equatorial plane of 
the Kerr-Taub-NUT spacetime j^] where circular geodesies do not exist. The remarks 
in 1^ about this spacetime ignored the equations of motion for 9 and are therefore 
incorrect when the NUT parameter £ is nonzero. 

The coordinate periods T± = ±27r/C(goo)± ~ 27r/|C(geo)±| of the two circular 
geodesies to complete one circuit with respect to the same static observer time 
line (with corresponding proper periods t± = j±^NT± as measured by the 
geodesies themselves) can be averaged to define an average coordinate period and 
its corresponding angular velocity 

^(0) = '^(C(gco)+ - C(gco)_) ' ^(0) = 27r/r(o) = 2/(C(g;|„)^ - C(gco)_) • (3-8) 
while an asymmetry parameter can be introduced by 

a = - (cr^ , +cr\ ) , (3-9) 

2^^(gco)+ ^(gco)_^ ' ' 

enabling one to express the C(gco)± in the form 

C(gco)± = l/(a±t^(o;) • (3.10) 

In the Kerr spacetime T(o) and lu(^q^ reduce to the Keplerian period and angular 
velocity, while the asymmetry parameter a reduces to the usual rotation parameter a; 
moreover a;(o) = {M/r^Y^'^, where M is the usual mass parameter and |a;(o)| = |ti| 
is the only nonvanishing Frenet-Serret torsion of these geodesies The Kerr ZAMO 
relative velocities and gamma factors are explicitly 

_ (1 T2aa;(o) )r^+a^ A^/^{1 ± au;^o^ ) 

AV2(«^(0)±1) ^(0)' '^±-i?(l-3A^/r±2ac.(o))i/2 ' ^ ^ 

where A = — 2A4r + is a factor appearing in the lapse function N = A^/^/i? 
and in the curvature ratio 

R/p = N{l-a^uj^^Q-^) (for Kerr only) . (3.12) 

With these definitions one can introduce naturally two angular velocity 
parameters and two time parameters, each of which has a physical interpretation. 
Table 1 summarizes these definitions. 

The average of the geodesic angular velocities defines the angular velocity of the 
GMPOs 

C(ginp) = (C(goo)+ +C(gco)-)/2 = -r''0t/r'"00 = ~gt4,,rl 9^4,,r 

= - au;^o)/(l - a^uof^^) < , (for a > 0) (3.13) 
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corresponding to a counter-rotation (if < a < 1) with respect to the static KiUing 
observers following the time coordinate lines. This angular velocity describes the GMP 
trajectories at all radii, regardless of the causal nature of the circular geodesies. The 
corresponding ZAMO relative velocity is the average of the geodesic relative velocities 

i^(gmp) = i^(C(gmp)) = 5(i^+ + i^-) • (3.14) 

Eq. (9.4) of 1^ shows that the GMP trajectories are 'Fermi- Walker relatively 
straight' with respect to the ZAMOs, which means that the direction of their relative 
velocity ±e^ in the ZAMO rest space appears to remain constant along their trajectory 
compared to the directions of nearby observers: P{n')^u{Cfi,^p^)^^ — 0, where P{n) is 
the projection into the ZAMO local rest space tangent to the time hypersurface. This 
in turn shows that at the GMPO horizon, the intrinsic Lie relative curvature |k((/>, nY\ 
of their circular orbit coincides with the single nonvanishing component \K{nY^\ of 
the extrinsic curvature of the time slice (the sign-reversed expansion tensor of the 
ZAMOs) along the orbit direction. 

By definition the GMP trajectories are respectively timelike, null or spacelike 

) takes the values 1, 0, —1. Let r-(g,jip) mark the radius 
at which they are null: the GMPO horizon. In the Kerr spacetime assuming a > 0, 
since C(gmp) < 0, these trajectories must become null before the time coordinate lines 
do as one approaches the horizon in the equatorial plane, so their observer horizon 
must occur outside the ergosphere (the time lines become null at the outer ergosphere 
radius) . 

Letting a tilde denote quantities defined for a complete geodesic circuit rather 
than a static observer circuit, since 

C(geo)± = C(gco)± - C(gmp) = ^2^C(gco)+ " C(goo)-) = ±'^(0)/(l " O^^fo)) , (3.15) 

the two geodesies have equal but opposite angular velocities C(goo)+ — ~C(goo)- 
with respect to the GMP trajectories. Thus if two distinct circular geodesies 
(C(geo)-i- 7^ C(goo)-) start at the same point on such a trajectory, they will return 
to it simultaneously, justifying the geodesic meeting point terminology. One such loop 
requires a coordinate time 

T{o) = 27r/[^(C(goo)+ - C(goo)-)] = ?'(o)(l - a^^fo)) , (3.16) 
which corresponds to a proper time orbital period 

■?(0) = 7(C(gmp))"^^T'(0) 

^\l-2M/r-a'^ujfft^{2 + M/r)\^^^Ti^o) (for Kerr only) (3.17) 
measured by the GMPOs and to 

(l-3M/r±2auj(o)Y/^ ~ , ^ , . 

= ^ ^^^(0) for Kerr only 3.18 

as measured by the geodesies themselves. The difference of the latter two proper 
periods is the observer-independent two-clock clock effect 

r+-f- = - lZ')Nf^o) = - 7l')7(C(gmp))r(0) ■ (3.19) 



when the sign e(C(gmp))^ of (2.7 
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null 


geos 


ZAMOs 


GMPOs 


c 


c± 

±1 


C(gco)± 
^(goo)± 


C(nmp) = 5(C+ + C-) 
i^(nmp) = 


''(gmp) ~ 2 (''(gco)+ ^" Qgco)-) 
'^(gmp) = +^-) 



Table 2. Tabic of angular velocities and corresponding ZAMO relative velocities. 



It is worth noting the following period ratio 

^(0) ^ ^ C(geo) + C(gco)- ^ ^ _ ^2^2 
T{o) (C(goo)+ - C(gco)-)^ 

R/p 



N 



(for Kerr only) (3.20) 



where the last equality follows from Eq. ( 3.12 ). 
Finally the product 

C(gmp)r(o) 271 ^'-^°°''^ ^ ^(gco)- ^ -2nauj(a) , (3.21) 

S,(gco)+ Mgco)_ 

is the change in coordinate angle cj) by which the GMPO has moved during one such 
circuit, itself the average of the individual angles that the geodesies have moved during 
this circuit 

C(gco)±T'(o) = 47r- = ±2^(1 T a^(o)) , (3.22) 

Mgco)+ Mgco)_ 

so aa;(o) and 1 =F aix>(o) are the number of circuits with respect to the static observers 
made by the GMPOs and geodesies respectively corresponding to each geodesic loop. 
For the Kerr spacetime, one may express the circuit number ati;(Q) of the GMPOs 
themselves in terms of the angular velocity ri(h) = a/(f+ + a^) of the black hole pO| 

auj^o)\r=r+ =A^r!(h)[2/(l-ar!(h))]i/' , (3.23) 

which is a proper fraction for < a/Al < 1 but has the value 1 for the extreme case 
a/Ai = 1 where the corotating circular geodesic becomes null at the horizon. 

Along a general circular orbit of angular velocity C starting at t = 0, = and 
parametrized by t, one has (j) = C^- The continuous parameter along the orbit is a 
parameter taking all real values in order to describe any number of circuits around 
the <?!>-coordinate circle in either direction. Its value mod 2tt gives a unique coordinate 
angle in the interval [0, 27r) for a point on such a curve. If one introduces a new angular 
coordinate (p = cj) — C(gmp)i which is comoving with respect to the GMPOs, then this 
induces a new parametrization (j) — {C ~ C(gmp))^ — C^i with : — > ±27r describing 
one complete circuit corotating or counterrotating with respect to the GMPOs. The 
two angular parametrizations are related to each other by 

^=(1-C(gmp)/C)'/'- (3.24) 

For the circular geodesies this becomes 

^ = (1 - C(gmp)/C(gco)±)0 = (1 T aw(o))(/) . (3.25) 

Table 2 summarizes the basic angular velocities and their corresponding ZAMO 
relative velocities instroduced so far. The average of the reciprocals of the geodesic 
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angular velocities defines the reciprocal angular velocity of the orthogonal direction 
within the t-cf) cylinder to another important observer family, the Carter observers 

1/C(car) = (C(geo)+ + C(geo)-)/2 = -T^^Pt/T^t = ~gt^,r/gu,r = a . (3.26) 

Note that in the static case gtc/i — 0, one has C(gmp) = = C(cai) = C(^ar)' -'^^ Kerr 
spacetime, the Carter observers and the ZAMOs share the same observer horizon, 
which is the (outer) horizon of the black hole. 

Another interesting observer family, the extremely accelerated observers p^ , 
extremizes the magnitude of the acceleration among the family of circularly rotating 
observers at each radius. Their angular velocity C(cxt) is given by a coordinate gamma 
factor weighted average of the geodesic angular velocities 

C(ext) = (r+C(geo)+ + r_C(geo)-)/(r+ + r_) . (3.27) 



4. Angular velocity maps 



The quadratic forms ( p.7|) and (3.5) may be associated with inner products on the 
space of spacetime angular velocity vectors dt + (d^ describing Killing trajectories. In 
both cases the condition of orthogonality with respect to the inner product 



(1 C2)(^^ ;^ j j =ACiC2 + s(Ci + C2) + c = o (4.1) 

leads to a fractional linear transformation between the two angular velocities 

B + ACi ACi ^ ' 

which is involutive: Ci ^ C2 — > Ci- 

For metric orthogonality {A,B,C) — (300, .90t, 5tt) this defines the bar map 

T C - C(th) g^o gtt 1 , . 

C =C(nmp)7 J — > 7 , (,4.3) 

<, - Xnmp) 5</'</' S 

giving the angular velocity of the direction orthogonal to the starting direction, where 
C(th) = —^/{g<j>t/gtt) is the angular velocity of the direction orthogonal to the static 
Killing observers. The bar map was already shown to be the reciprocal map when 
expressed in terms of ZAMO relative velocity in (|]^). Note that the null directions 
are invariant under this reciprocal map, so <^± = C±- 

For vanishing radial covariant derivative of one with respect to the other, 

- 2[Vo,+co, [dt + Z{Od^)]r = AC Z(C) + 2B[C + Z(C)] + C = , (4.4) 



which is the quadratic form corresponding to the condition (^^) satisfied by the 
geodesic angular velocities, then {A,B,C) — {g(t>(t>,r,g<i>t,r,gu,r) and this defines a new 
map Z 

/■ C - C(car) <? j,t-^ gtt.r 1 , . 

2^(0 = C(gmp)7— 7 > T, (4.5) 

S ^(gmp) y(f>(t>,r S, 

satisfying Z{Z{C)) — (. Z{C) is the angular velocity of the stationary axisymmetric 
vector field tangent to the symmetry orbit which is covariant constant along the Killing 
trajectory with angular velocity C. 
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Evaluating the map Z in terms of the ZAMO relative velocity, one finds 

.„,^^^ '^{0 - >^Z{nmp) 9,#.t— 

0) = i^(gmp)— TTT ^ ' -I'+V- VKQ . 

>^{0 - '^(gmp) 

7(gmp)k ~ '^(ginp)l „\ 

= |[.~.(Z(C,))][.-.(^(C-))]|V^ ' ^'-'^ 

where the vanishing ZAMO relative velocity i^(C(nmp)) = of the ZAMOs themselves 
shows that 

i^Z(nmp) = l^(Z(C(„n,p))) = + yZ^Vl2 , (4.7) 

the explicit value following from the actual calculation. 

Since Z{Z{C,)) = C, the vector with angular velocity Z{C,) is null at the angular 
velocities 

2(C±)=C(gmp);^— (4.8) 
S± S(gmp) 

or equivalently in relative velocity form 

^^l-ZlCi) = -; ; = -z — > ^v+v- . (4.9) 

From this one easily derives the inequalities 

< v{Z{C-)) < v{Z{C+)) < v+ , (4.10) 

with equality holding for the outer inequalities in the separate limits i>± —>■ ±1 when 
the corresponding geodesies become null. These outer inequalities separately reverse 
when the corresponding geodesies are spacelike. Thus in the interval Z{(-) < C < 
Z(^+) or equivalently iy{Z{(-)) < < v{Z{(+)), the vector dt + Z{()d^ which is 
parallel transported along the curve with angular velocity ^ is spacelike, and outside 
this interval it is null (at the endpoints) and then timelike. The bar map sends this 



vector to an orthogonal one whose 1-form according to (2.10) is proportional to 



YiZiOf ^R[dcj,~Z{Odt] , (4.11) 

which is then timelike within this interval and spacelike outside it, while remaining 
null at the endpoints. 

5. Closed circular orbit holonomy 

Parallel transport of a vector X around the (/)-coordinate circles (parametrized by 
the coordinate 0) is governed by the constant coefficient linear system of differential 
equations for the coordinate components 

d X" 

V^X" = 0^— — -A%X^ = 0, (5.1) 

which has the general solution 

X(0) = e*'^X(0) . (5.2) 

The matrix = — r"^^ = —g°'^gct,[f3,fj,], which consists of the coordinate components 
of a mixed tensor A whose covariant components Aap — —Ap^ are antisymmetric and 
therefore define a 2-form, generates a 1-parameter family of Lorentz transformations 
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of the tangent space expressed in nonorthonornial coordinates on that space. This 
matrix (constant along the ^-coordinate circles) evaluates explicitly to 

A% = (lni?),,[dr A r(C(gmp))']% = R/P [4 A r(C(g,np))1% . (5.3) 

Since this is the mixed index form of a decomposable 2-form (wedge product of two 
1-forms), interpretating its exponential as a Lorentz transformation is simple. When 
nonnull, decomposing this 2-form into its magnitude a and a unit 2-form B (with 
unit forms suitably defined) yields the orientation of the 2-plane (spanned by the 
corresponding vectors) in which either a boost (timelike case) or rotation (spacelike 
case) takes place and the corresponding pseudoangle or angle is that magnitude a, 
while in the null case, normalization is not possible and a null rotation takes place in 
this plane. The following factorization does just that 

A% = = (7(0)^% , (5.4) 

where recalling the notation of ( ^.10[ ) one defines 

O- = 7(C(gmp))"^-R/P , cr(0) = R/P , 

B'^e'.A i7(C(gmp))' , B'^e^A y(C(gmp))' , (5.5) 
with the properties 

= -e(C(gmp))'S , S3--e(C(gmp))'7(C(gmp))-'S. (5.6) 

Recall that the unit vector C/(C(gmp)) when spacelike is just the (f) direction in the 
GMPO local rest space. 

Thus the 2-plane of this Lorentz transformation is spanned by the radial direction 
and the direction within the t-cf) cylinder orthogonal to the GMP trajectory direction, 
with the factor e(C(g,jip))^ — 1, 0, —1 describing these three cases in which the 2-plane 
is spacelike, null, and timelike respectively, corresponding to a rotation, null rotation, 
and boost. Thus outside the GMP horizon in a Kerr spacetime, this is a rotation in 
the r-(/) plane in the GMPO local rest space and the identity in the orthogonal 2-plane 
spanned by [^(C(gmp)) and eg. 

The ex pon entiation is easily evaluated by collapsing the infinite series using the 
properties ( ^.6[ ); letting 1 stand for the identity matrix in the exponentiation, one 
finds 

r^'P'TB ^ [l-7']+COs[e(C(g„,p))a0]P 

e"^^ = J +e(C(gmp))"^ sin[e(C(gmp))o-0] B , e(C(gmp)) ^ , (5 7) 

[ e^^^"'" = l + a^o)4>B+laf,)<l,'B^ , e(C(gmp)) = 0, 

where the matrix / mixed tensor 

r - -e(C(gmp))"'i?' , e(C(gmp)) ^ , (5.8) 

satisfying — V and BV = B — VB, is the projection into the 2-plane spanned by ef 
and U (C(gmp) ) • The exponential matrix represents a Lorentz transformation acting in 
this 2-plane in the nonnull case, while acting as the identity in the orthogonal 2-plane 
with projector l — V. When e(C(gmp)) = 1) it is a clockwise (negative) rotation by the 
angle acj) with respect to the ordered orthonormal axes (e^ , C/(C(gmp))) of the boosted 
spherical frame. Note that dt + C(gmp)f^</' is the direction of the 1-parameter family of 
stationary axisymmetric vectors tangent to the t-4> cylinder which is invariant under 
this transformation and hence is covariant constant along the curve. 
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The eigenvalues of the matrix A are 0,0, ±/, where the frequency function 
/ = £(C(gmp))<'' can be real, zero or purely imaginary. One may also obtain decoupled 
equations by differentiation and backsubstitution instead of the eigenvalue-eigenvector 
approach (equivalent to the above explicit exponentiation). Thus 

= a^[B^xr = -f^iaBXf = -/'^ (5.9) 



leads to decoupled equations for each of the der ivatives dX"/d(j), for example. For the 
Kerr spacetime, setting a = t leads to Eq. (LI) used by Maartens et al ||]. 



Note that for asymptotically flat spacetimes like Kerr, as r ^ cxo one has the limits 
{grr, gci><t>) (l;''^)j P ^ 7(C(gmp)) ^ 1, and (T 1 leading to an angular velocity 
of 1 in the r-(j) plane of the tangent space, while B generates a clockwise rotation of 
the ordered boosted r and (j) axes. This causes a parallel transported vector to rotate 
in this plane by an angle in the clockwise direction with respect to the boosted 
coordinate axes to maintain its direction as those axes rotate in the counterclockwise 
direction by the same angle as (/) increases. One must subtract this orbital rotation by 
the angle <j} in order to obtain the rotation of the parallel transported axes with respect 
to the nonrotating axes at spatial infinity, corresponding to the angular velocity cr — 1 , 
or 1 — (7 in the usual counterclockwise direction. 

The dimensionless frequency function can be expressed in the form (see Eq. ( 3.20D ) 



~ R/p ^ ( R/p ^(0) 



1 2X 2 2 

1- — -«H)(2 + 



1/2 

(for Kerr only) . (5.10) 



It converts the angle of revolution of the geodesic with respect to the static observers 
to the corresponding angle that a parallel transported vector rotates in the opposite 
direction with respect to the spherical frame about a fixed axis in the GMPO local 
rest space, and in the Kerr spacetime it is just the ratio of the GMP proper geodesic 
period to the average static observer coordinate period for the geodesic pair. 
In the combination 

= 7(C(gmp))"M^0]/P , (5.11) 

representing the total parallel transport angle from to ^ (relative to the spherical 
frame), the bracketed expression is the arclength along the coordinate circle and 
the gamma factor converts this to the arclength along the corresponding circular 
trajectory orthogonal to the GMP trajectories (Lorentz contraction when timelike) 
which projects down to the coordinate circle along the ZAMO trajectories, and finally 
the last factor takes its ratio with respect to the local radius of curvature of the 
coordinate circle as seen by the ZAMOs. The angular velocity goes to zero exactly 
where the gamma factor goes to infinity at r = r(gi„p), the horizon for the GMPOs. 
This relation can be restated in terms of the embedding deficit angle as 

^ = 7(C(gmp))-^(l-%^) . (5.12) 
Then for a change of angle <j) — ±27r corresponding to one revolution, one has 

= ± 7(C(gmp))"^(27r - A(dof)) , 
(a - 1)0 = ± 27r[7(C(g„,p))"i - 1] T 7(C(gmp))-' A(def) . (5.13) 
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Since B generates a clockwise (negative) rotation with respect to the ordered axes 
(cf , f7(C(gmp))) of the boosted spherical frame, the corotating/counterrotating circuit 
leads to a counterclockwise/clockwise rotation of a paraUel transported vector with 
respect to these two axes by the deficit angle Ajjcf) when the gamma factor is close 
to 1, corresponding to Figure 2 of Thorne 

Where both circular geodesies are timelike and a is positive and has a rational 
value a = p/q at a given radius, where p and q are positive integers, then after q 
revolutions ((/> — i:2T:q), any vector undergoes a rotation by a4> = ±27rp and will 
thus rotate back to its original starting value ('holonomy invariance') so the holonomy 
group restricted to this family of curves (fixed p, variable q) must reduce to a discrete 
subgroup of the rotation group. For the Kerr spacetime ct is a monotonically increasing 
function of r from its zero value at the GMPO horizon ?'(gmp) to its asymptotic limiting 
value 1 at r — > oo, so if cr = p/q, then p < q and q > I. This means that a transported 
vector undergoes fewer revolutions with respect to the boosted spherical frame than 
the number of loops around the hole before returning to its original orientation. The 
radii at which such holonomy invariance occurs for the 'harmonic' parameter family 
q = 2,3, . . ., p = 1,2, . . . ,q ~ 1 cover the interval (r(gmp), oo) densely in the same way 
that the rational proper fractions cover the interval (0, 1), creating the band structure 
referred to in ||]. 



6. General circular orbit transport 

In order to compare the parallel transport of a vector after it returns to the same 
GMP trajectory after being carried along by the two oppositely-rotating circular 
geodesies, one first needs to extend this discussion to circular orbits with constant 
angular velocity ( £ R. These curves can be parametrized naturally by either the time 
coordinate t or when C 7^ by the extended values (j) = (^t oi the (/^-coordinate, assuming 
such a parametrized curve starts at t = and (ji — for simplicity. Corotating and 
counterrotating orbits (of any number of circuits with respect to the reference orbits 
of the static observers with C = following the time coordinate lines) correspond to 
positive and negative angular velocities respectively. Note that the previous case of 0- 
coordinate circles corresponds to the hmits lim^_»±oo C~^{9t + (d^) — d^. One can also 
use a continuous parametrization corresponding to the comoving angular coordinate 
^ = (1 — C(gmp)/C)0 relative to the GMPOs as a parameter along this orbit. 

The equation Vt/(n^ = of parallel transport along such a parametrized curve 
has the same form as ( |5.l| ) 

{d 
-— - A{CT,,X^^ = , C 7^ , 

,?a (6-1) 
QA(0\X^ = , C = , 

where the matrix 

A{0\ = -(r%M + C-'r",^) = -C'9"^\gt[M + (:9m^^]) (6-2) 

is again antisymmetric when index-lowered and so corresponds to the mixed tensor 
form of a 2-form and generates a Lorentz transformation, although expressed in 
nonorthonormal coordinates on the tangent space. It is obtained from A in the new 
derivation adding an extra term to ^(C(gmp))'' = R['^4> ~ C(gmp)di] in ( pT^ ) yielding 
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instead 

C(gmp)dt - (C(gmp)/C)(d'/> - C(car)dt)] 

gmp) ; 

C(gmp))-Rdt , C = C(gmp) , (6-3) 

ap)/C(car))^d(^ , C = C(car) , 

thus interpolating between (whe n C = C(car)) and dt (when ( = C(gmp))- 

The angular velocity map Z of ( |4.5| ) determines the trajectories orthogonal to the 
plane of the Lorentz transformation. Direct calculation shows that it satisfies 

^(C(gco)±) = C(goo)± , 2^(C(oxt)) = C(oxt) , ^(C(cxt)) = C(oxt) , 

2{0) = C(car) , 2(C(car)) = = Z(C(g„,p))-l , ^}^f^ZiO = C(gmp) , (6.4) 



where C(ext) is the angular velocity of the extremely accelerated observers (B.27). The 
limits C ±oo reduce the present case to the previous case for the (/)-coordinate circles, 
while for C = this shows that along the time coordinate lines, the Carter observer 
local angular trajectories play the same role that the GMP trajectories play for the 
(ji-coordinate circles. The boost from U{() to U{Z{C,)) when both are timelike maps 
the local rest space of the orbit along U{Q to the local rest space of the plane in which 
the parallel transport Lorentz transformation takes place. The vector U{Z{()) is a 
unit vector along the (/)-direction in this latter local rest space and {ef,eg,U{Z{())} 
is the boosted spherical frame. Figure 1 shows the orientation of the intersection of 
the Lorentz transformation plane with the t-(f> tangent plane (i.e. U{Z{Q) — U{Z{Q) 
modulo sign) for the directions of the various 4- velocities which arise in this discussion. 
The exponentiation for the solution of the parallel transport equations 

X{c^) = e^^^^^X{Q) (6.5) 

and its interpretation as a Lorentz transformation in a 2-plane is again accomplished 
with the relations 

A(C) = a(C)i?(C)=a(o)(C)S(C) , 

<y{0 = l{Z{0)-\R/p] (1 - C(gmp)/C) , fT(o)(C) = [R/p] (1 - C(gmp)/C) , 

B{0' = A U{Z{0)' , BiCf = 4 A ^(^(C))' , 

BiCf = -e{Z{OfB{0 , = -6(Z(C))^(Z(C))-2S(C) . (6.6) 

The general exponential solution then has exactly the same form as (5.7), in terms 
of the ^-generalized quantities 

0'x(C)B(c) ^ [1 _ ^ cos[6(2(C))fT(C)0] V{0 

e^^(C) = <| +e(Z(C))-isinKZ(C))a(C)0]B(C), e{Z{0)^0, (g^) 

g0.,o,(C)B(C) = 1 + a(o)(C)0S(C) + \<T(,){Cfcf>''B{(:f , 6(Z(C)) =0, 
where now 

ViC) = ~e{Z{C))-'B{Cr , e{Z{0) ^ . (6.8) 

Note that by defining the new dimensionless frequency with respect to the GMPOs 
by (t(C)0 = (t(C)(/', one finds 

a(C) -(i-C(gmp)/C)-V(C)=7(2(C))-'i?/p, ^(2(0) 7^0, 

^(o)(C) = (1 - C(gmp)/C)"'^(o)(C) = RIp , e(2(C)) = , (6.9) 
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Figure 1. The orientation of the intersection of the parallel transport Lorentz 
transformation plane with the t-<j) tangent plane for the various directions inside 
the half circle for a Kerr black hole with a/J\4 = .5 at r/M = 4. The horizontal 
and vertical directions are along and n (ZAMO angular direction and 4- 
vclocity) respectively, while the null directions (dashed lines) are at 45 degrees 
from the horizontal. Each arrow inside the half circle indicates a direction in the 
t-rf) velocity space and the direction arrow attached to it indicates the direction 
given by the composed map Z applied to the corresponding relative velocity. 
The direction of the covariant constant vectors for a given initial direction is the 
direction orthogonal to the attached arrow in the Lorentz geometry (reflection 
across the 45 degree line). The double arc indicates the boost zone. The extremely 
accelerated observer direction and its orthogonal counterpart indicated by ext 
and ext are the only ones aligned with their boost /rotation plane direction, while 
the two geodesies are the only ones aligned with their covariant constant vector 
direction. The maps Z and Z are actually both reflection symmetric across the 
extremely accelerated observer direction. 



SO it is the angle with respect to the GMPOs which plays the same role as the 
coordinate angle for the 0-coordinate circle discussion. When Z{C,) corresponds to a 
timelike direction, the inverse gamma factor in a locally Lorentz contracts spacetime 
arclength along the 0-coordinate circle to spacetime arclength along the corresponding 
tilted local rest space direction. 

The local direction of the plane of the Lorentz transformation within the t-4> 
cylinder is now along C7(Z(C)) (when e(Z(C)) ^ 0) or Y{Z{C)) (when e(Z(C)) = 0). 
For the timelike circular geodesies, which arc invariant under the map Z, this is just 
the angular direction in their local rest space, while along the time lines this is the 
local time direction of the Carter observers. In general when the Z{() trajectories are 
timelike/spacelike, this represents a rotation/boost; the horizon for the corresponding 
observers is reached when 'y{Z{())~^ vanishes. This occurs at the two angular 
velocities C = -2(C±) satisfying the inequalities (4.5). At a given radius, the interior of 
this angular velocity interval represents the boost zone, while its exterior represents 
the rotation zone, separated by null rotations which occur at the endpoints. 

Recall that dt + Z{C,)dtjj is the direction of the stationary axisymmetric vectors 
which are tangent to the t-(j) cylinder and invariant under this Lorentz transformation 
and hence are covariant constant along the curve. The extremely accelerated observer 
world lines are the only circular orbits for which this covariant constant angular 
direction is orthogonal to the orbits themselves, corresponding to the phase-locking of 
Fermi- Walker transported spin vectors found by de Felice and Usseglio-Tomasset |l2| . 
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In fact since Z must be a fractional linear transformation of the relative velocity with 
respect to any observer (since the velocity addition boost formula is a fractional linear 



transformation), the four conditions on the first line of (6.4) show that with respect 
to the extremely accelerated observers 2 : v —v^v^/v (where = — i^+). Thus 
these observers see both the geodesic relative velocities and the map Z symmetrically 
in the relative speed 

If one considers circular orbits with a constant (finite) angular velocity, the orbits 
are not closed and one cannot consider spacetime holonomy. One can obtain a closed 
loop by joining together two distinct orbits of this type at a crossing point, which 
then periodically meet each other, forming a sequence of closed loops in spacetime. 
For example, in the Kerr spacetime one could choose the first orbit to have an arbitrary 
angular velocity ^ and the second orbit to be a time coordinate line with zero angular 
velocity, which is a world line of a static observer, along which vectors are boosted 
by parallel transport. The net holonomy transformation would then be the inverse of 
one Lorentz transformation times the other Lorentz transformation. Rothman et al 
slide over this point by moving quickly on to the case in which the two orbits are 
oppositely-rotating geodesies. 

It should be noted that Marck has explicitly constructed a parallel transported 
frame for general geodesies in the Kerr black hole using its Killing- Yano tensor, a 
result which extends to a wider class of stationary axisymmetric spacetimes |l^ . 

7. Timelike geodesic holonomy 

Along the timelike geodesies starting at t = and cj) = and parametrized by t, 
the parametrizations by the GMPO comoving angle and the coordinate angle and the 
corresponding dimcnsionlcss frequencies are now related to each other by 

= (1 - C(gmp)/C(gco)±) = (1 T aW(o))"^ = ^(C(geo)±) ^ 
9± C^lS(gco)±j 

while 

HC(gco)±) = i±'R/p ■ (7.2) 

It is convenient to introduce the abbreviations 

0'± = 0-(C(goo)±) , O-(avg) = (0'++0--)/2 , (T(dif) = ((T+-(T_) , /| = e(C(gco) ± )^CT± (7.3) 

and their corresponding tilde versions. The quantities are the frequency-squared 
functions appearing in the decoupled second-order equations analogous to ( ^.9[ ) for 
parallel transport along the circular geodesies. For the Kerr spacetime because of the 



special condition (3.20), one has 



II 37Wio |i/2 ^± 

^± = |1- — ±2au;(o)|'/'(lTa^(o)) = 7^ ■ (7.4) 

Thus in the Kerr spacetime the dimensionless frequency factors which convert the 
angle of revolution with respect to the GMPOs to the parallel transport angle with 
respect to the boosted spherical frame are just the ratios of the individual geodesic 
proper periods to the average coordinate period with respect to the static observers. 
The corrected frequency factor a± — l corresponds to the rotational angle with respect 
to boosted spherical axes whose orbital rotation is removed & p|, O . 
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Figure 2. Radius ratio r/M versus a/M in the Kerr equatorial plane for 
the horizon radii of the ZAMOs (hor), the static observers (erg), the GMPOs 
(gmp), and the counterrotating and corotating geodesic observers (geo— , geo+). 
As a/A4 increases from to 1, the horizon radius ratio r/J\4 decreases from 2 
to 1, the ergosphere ratio remains at 2, and the GMP horizon ratio increases 
from 2 to 1 + V2. The curves marked 'geo—' and 'geo+' are the radii of the 
counterrotating and corotating null circular geodesies, which start at 3 and end 
at 4 and 1 respectively. 

These frequency functions {a± and a±) go to zero when the corresponding 
geodesies become nuU as 0, corresponding to the single GMP horizon in 

the case of the corresponding discussion for the 0-coordinate circles. For the Kerr 
spacetime, these 'geodesic horizons' mark the regions in the equatorial plane where as 
one approaches the black hole, first the counterrotating geodesic goes null (at r = rAB, 
where cr_ = 0) and then the corotating geodesic goes null, (at r — rsc, where cr+ = 0), 
which divides the equatorial plane into three regions A, B, C, where of the two circular 
geodesies respectively both are timelike or one or neither is timelike. Each such radius 
represents the only zero of the corresponding frequency function outside the horizon. 
Explicit formulas for these two radii are given by Chandrasekhar Figure 4 of |^ 
shows a plot of these radii and the horizon and ergosphere radii as a function of a for 
the physical interval a G [0,1]. Figure 2 reproduces these curves together with the 
GMPO horizon radius ?'(gmp) = f*- 

When both circular geodesies are timelike, the parallel transport Lorentz 
transformation is a rotation by the respective angles o'±(/)±, but this rotation takes 
place in different planes, namely in the local rest space of each geodesic. The geodesies 
meet at every half-revolution with respect to the GMPOs; the net parallel transport 
angles along each are 

zbTTCTi = ±7r7±i[i?/p] . (7.5) 

Apart from the gamma factor, this corresponds exactly to the shortened rotational 
angle by which a transported vector rotates in the flat tangent cone to the embedding 
surface during each half revolution with respect to the ZAMOs, and the gamma factor 
then Lorentz contracts this amount to correspond to the angle in the geodesic local 
rest space which projects down onto this angle in the ZAMO local rest space. 
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After a positive integer q of such half-revolutions, at either the original GMPO 
world line {q even) or the antipodal GMPO world line {q odd), the net parallel 
transport angles along each are 

±qTTa± = ±qTT^^^[R/ p] . (7-6) 

Joining two oppositely-rotating geodesic curves together at their common departure 
point (initial meeting point), one obtains a closed loop representing q complete circuits 
with respect to the GMPOs. With a common initial vector X{Q) at this starting point 
at = = (/), the two parallel transported vectors X{(j)±) at the final meeting point 
will differ by a holonomy transformation belonging to the Lorentz group. 

The orthogonal projections of X((/)) into the two local rest spaces of the pair of 
oppositely-rotating geodesies define its apparent directions as seen by those geodesies. 
Only this orthogonal component evolves along each geodesic, undergoing a rotation 
in a fixed (coordinate component) plane. The holonomy transformation at the final 
meeting point then results from applying the inverse of one of these rotations, followed 
by the boost from that local rest space to the other, and then applying the other 
rotation. At those radii for which this transformation is the identity matrix, one finds 
holonomy invariance. Holonomy invariance for the entire tangent space requires that 
each rotation separately return to the identity, since the two rotations take place in 
different planes and the only way one can reassemble the two pairs of orthogonal pieces 
of the final vectors X((j)±) after parallel transport to form the same spacetime vector 
if the temporal parts are unchanged is if the spatial parts are also unchanged. 

In the static case, where C(gmp) = and a± = ij±, the geodesic gamma factors 
7+ = 7_ are the same and the two rotations in opposite senses by the same angle in 
absolute value will agree if they are each equal to a rotation by some multiple of 27r, 
so 

7±^[i?/p] gTT = 27rp or a± = jg'^R/p = — , (7.7) 

where p is some positive integer representing the number of parallel transport loops 
which take the the initial vector to the final vector along each such curve separately 
(with respect to the boosted spherical frame along the way). This occurs when 
the GMP frequency functions a± assume rational values of the form Ipjq. If 
Rj p < 1 as in the static Kerr case (Schwarzschild), one must have 2p < q and 
q > 2. Here the formula (|3.1l|) reduces to 7± = (1 - 2M/r)'^/'^/{l - 3X/r)i/^ 
thus converting the 7(C(gmp)~^-R/p = R/ P = (1 ^ 2M/r)^/-^ factor for i/i-coordinate 
orbits to ^^^R/p = (1 — SM/rY^^ for geodesies, recovering both results of Rothman 
et al |l[] for the frequency functions, but they neglect the crucial factor of 2 which 
rules out the odd values of their parameter m = 2p. Note that 2p is the total number 
of loops with respect to the boosted spherical frame that the final vector AT ((/>_) 
undergoes as it is parallel transported back from the final meeting point to A(0) 
along q counterrotating half-circuits and then on to the final vector A(0+) after q 
corotating half-circuits returning to the final meeting point. In this way the 'even' 
members of the harmonic parametrization q — l,2,...,p = 1,2,. ..,(7 — 1 generate a 
corresponding discrete spectrum of radii at which the holonomy subgroups restricted 
to these orbits are finite. These radii densely cover the radial interval where the two 
circular geodesies are timelike, reproducing the band structure of holonomy invariance 
found for the closed 0-coordinate loops. 
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These geodesic gamma factors are distinct in the nonstatic case (c7+ ^ c7_), 
whore the condition that the separate rotations each return to the identity leads to 
two independent conditions rather than one as in the static case 

a± qn = 2ttp± -> a± ^ 2p±/q , (7.8) 

where now p+ ^ p-. In general it is not possible to satisfy both conditions at the 
same radius, although exceptional solutions of this type can occur at a set of radii of 
measure zero in the interval where both geodesies are timelikc. This is clearly seen 
in the Kerr spacetime case (see Fig. 2) where as one increases the rotation parameter 
a away from 0, the graphs of a± versus r move away on opposite sides from the 
Schwarzschild graph. The counterrotating geodesic frequency curve moves to the 
right like the corresponding horizon rAB of the counterrotating geodesies, while the 
cororating geodesic frequency curve moves to the left like the corresponding horizon 
tbc of the corotating geodesies, with their zero values marking these two respective 
radii. For fixed q and a given pair p-^. < p_ in Schwarzschild, the corresponding radii 
''(9,p±) which a± = 1p±/q satisfy T"(g^p^) < r(g^p_), but as one increases the rotation 
parameter a for a given the various radii pairs eventually cross, one by one, modulo 
large a effects. However, for a given value of a, there is no guarantee that even one 
such common radius exists. Thus the band structure of holonomy invariance which 
exists for the static case does not exist in the nonstatic case, where it has instead 
been diluted by spreading out over the parameter interval of a values, unlike the case 
p+ = P-, which all occur at the single value a = 0. A similar situation probably arises 
for other spacetimes of this elock-efFcct family. 

Figure 3 illustrates this situation for the Kerr frequency functions with a/ = 0.5 
compared to the Schwarzschild case (middle curve). The horizontal lines at 3:, |, | 
illustrate either a q = 4 circuit loop and 2p± = 1,2, 3 or ag = 8 circuit loop 
and 2p± = 2,4,6. For example, the intersection points with the lines at 2p±/8 
mark the radii at which the parallel transport rotation separately returns to the 
identity compared to the spherical frame after an 8-cireuit loop, in the process making 
p± = 1,2,3 revolutions with respect to that frame. One can see that the pair 
(p+,p_) = (2, 1) has already crossed at a smaller value of o, while the pairs (3, 1) and 
(3, 2) have yet to cross for larger values of a, all of which are easily found numerically. 

The real problem here is the boost between the local rest spaces of the pair of 
geodesies. If one were interested in comparing Fermi- Walker transport along curves 
modulo such boosts, then an analog of holonomy invariance still exists. Fermi- Walker 
transport along a timelike curve describes the behavior of the spin of a test gyroscope 
carried along such a world line, and reduces to parallel transport along a geodesic. 
Spin vectors lie in the local rest space of the world line along which they travel, so to 
compare, at a spacetime point, spin vectors carried by observers in relative motion, 
one has to boost the spin vector from one observer to the other. 

The degree to which this is physically interesting can be disputed, but if one 
accepts it mathematically, then one can parallel transport a spin vector around a closed 
geodesic loop like a baton in a relay race, passing off the baton from one geodesic to 
the other at the start and finish by a boost in order to make the comparison in the 
orientation of the spin direction without the relative motion distortion. For the case 
at hand, the boost comparison identifies the angular direction along (j) in the local rest 
spaces (namely f/(C±)), allowing one to ignore the fact that the rotations take place 
in different planes. 
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Figure 3. The GMP frequency functions a-± versus r/JV4 in the Kerr equatorial 
plane for the corotating and counterrotating geodesic observers (geo+, geo— ) in 
the the Kerr spacetime with a/M = 0.5, together with the a = limit (middle 
curve). The black hole horizon is at r/A4 = 2. Intersections with the three equally 
spaced horizontal lines between and 1 mark the radii at which the parallel 
transport rotations separately return to the identity after q = 8 GMP circuits 
corresponding to p = 1,2,3 (left to right) complete revolutions with respect to 
the spherical frame. 



On the other hand one can imagine an experimental situation where each spin 
system is interrogated separately at the crossing event. Then the final rotations will 
agree (ignoring the boost) if the two separate angles are the same mod 2tt after g > 
separate half-circuits (joint whole circuits) 

P 

97rCT(gco)+ - (-g7ra-(goo)-) = 2g7r(7(avg) = 2np -> CT(avg) = - , (7.9) 

i.e., when the average frequency function o'(avg) assumes rational values. This leads 
to a situation analogous to the static case, but now including the odd values of the 
numerator p. The average of the two individual angles instead gives the relative angle 
of their common direction (modulo the boost) with respect to the spherical frame at 
the meeting point 

^[97rCT(geo)+ + (-'77rCT(geo)-)] = ^'77ra(dif) = ^g7rg ^+, , (7.10) 
Z ^ r(o) 

which is directly proportional to the fractional two-clock clock effect given by the ratio 
of ( [3.1S| ) and ( 3.17 ). The proportionality factor is just the frequency function a for 
the (/i-coordinate circle discussion multiplied by ^(/tt. 

If one imagines a gedanken dual drag-free satellite gyroscope experiment with 
two oppositely-rotating equatorial plane geodesic orbits at the same radius, then it 
makes sense to assume that the internal spin up mechanisms would give the gyros the 
same spin magnitude in their rest frame relative to the satellite housing itself, which is 
presumably locked onto a distant guide star and hence to the local spherical frame in 
a calculable way. Since each spin vector on these orbits precesses around an axis fixed 
with respect to the boosted spherical frame, one could compare the total angles of 
precession at meeting points with each other as well as with the local spherical frame. 
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in which case the above mathematical considerations would come into play, including 
the use of boosts to compare orientations of gyros and axes in relative motion. These 
individual rotations of the spin vectors represent the combined geodetic and Lense- 
Thirring precession effects, which are aligned for these orbits, in contrast with slowly 
precessing polar orbits in the limit of small a where they are orthogonal, as in the 
GP-B experiment . 

8. Conclusions 

Previous work on the holonomy associated with circular orbits in particular stationary 
axisymmetric spacetimes has been generalized to all those which admit clock effects, 
using the Kerr spacetime as an illustrative example. The details of the calculations 
depend heavily on the geometry of the GMPOs which are themselves defined by 
the geometry of clock effects. They directly determine the parallel transport 
transformations relative to the natural symmetry-adapted frame along the closed 
curves tangent to the Killing vector associated with axial symmetry, and do the same 
in combination with the Carter observers for general circular orbits. In the nonstatic 
case the band of holonomy invariance which exists for special geodesic clock-effect loops 
in the static case disappears. With some modifications, one can extend these parallel 
transport results to the case of spin precession to compare the total precession angles. 
In the Kerr spacetime the dimensionless rotational frequencies associated with parallel 
transport rotations are particularly simple ratios of various circular geodesic period 
quantities. The rich mathematical structure that has emerged from this geometric 
problem has been made apparent only through the relative observer analysis associated 
with gravitoelectromagnetism. 
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Appendix A. Prenet-Serret approach 



Accelerated curves have a natural orthonormal frame, the Frenet-Serret frame, with 
respect to which the parallel transport Lorentz transformation is described in terms 
of the curvature k (magnitude of acceleration) and the first and second torsions ti 
and T2 (Fermi- Walker rotation parameters). The relation between this analysis and 
the coordinate frame analysis performed here is rather simple. 



The general Frenet-Serret frame along a circular orbit 1 18, 3|, |19, 20 , for which 
Eq = U{Q, satisfies 



^ - F^aEp , = [nEo A El- n^i A E2 - T2E2 A Es^^ 
ar 



(A.l) 



which is explicitly 



DEn 



= kEi 



DE^ 



= kEq + T1E2 



DE, 



= -TiEi 



T2E3 



DE-, 



dr dr dr " dr 

A vector Z = Z°'Ea that is parallel transported along U then satisfies 
DZ" dZ 



= -T2E2 . (A.2) 



dr 



dr 



(A.3) 
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For this case the frame is given by Eq = C/(C) , Ei = Cf , E2 = U{() , E^ = ' ^^'^ 
the curvature and torsions are 

^ - -p-^(C)'(KC) - ^+)(KC) - 1^-) , n = -^7(C)-'§ , r2 = , (A.4) 

leading to the Frenet-Serret relations 



From the definitions (6.1) and (|6.q), which are equivalent to 



^ = aiOtJiZiO) , (A.6) 

and the second Frenet-Serret relation together with the reparametrization relation 
d(j)/dT = jiO^/N, one has, when U{C) and U{Z{C)) are nonnuU, 

C7(Z(C)) = ± |rf - + r,U{0] , 

C/(Z(C)) = ± |rf - kY'^'[t,U{0 + ^UiO] , 

aiChiCK/N = ± - k2|1/2 = ±|i^-/3^^„^/2| V2 ^ (a.7) 
modulo signs which depend on the signs of k, ti, C and where 

Ixl . \(u u ^1l/2 7(C)7(t/(^(C)),t/(C)) .... 
P 7(gmp) 



Modulo these signs the first two lines of ( A.7 ) represent a boost from (U{C)i U{C)) to 



(J7(Z(C)), U{Z{())) with relative velocity k/ti when \k/ti\ < 1, and a boost plus the 
bar map when | k/ti | > 1 . 

The case of Fermi- Walker transport instead of parallel transport correspon ds to 
the replacements D £'(fw)j k — > in the above relations, while projection of ( |A.6|) 
into the local rest space of U{C) yields the Fermi- Walker derivative along U{C) of the 
vector field Cf , which is spatial with respect to all circularly rotating observers 

= P(UiO)^ = ^{OP{U{C))U{Z{0) = a(C)7(C/(2(C)), UiO) UiO . (A.9) 
Thus the Fermi- Walker transport dimensionless frequency is instead 

'^(fw)(C)=7(f/(2(C)),f/(C))^(C) , (A.io) 

with an extra relative gamma factor describing the inverse Lorentz contraction which 
occurs in the projection, which has unit value for geodesies where this projection 
reduces to the identity. The relative gamma factor may be expressed in terms of the 
ZAMO gamma factors |l| as 7(C/(2(C)), U{C)) = 7(2(C))7(C)[1 - HZ{C)HC)]. 
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